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ABSTRACT 

The dynamics of vibrations in flexible structures can be conveniently modeled in terms of 
frequency response models. For structural control such models capture the distributed param- 
eter dynamics of the elastic structural response as an irrational transfer function. For most 
flexible structures arising in aerospace applications the irrational transfer functions which 
arise are of a special class of pseudo-meromorphic functions which have only a finite num- 
ber of right half plane poles. In this paper, we demonstrate computational algorithms for 
design of multiloop control laws for such models based on optimal Wiener-Hopf control of 
the frequency responses. The algorithms employ a sampled-data representation of irrational 
transfer functions which is novel and particularly attractive for numerical computation. One 
key algorithm for the solution of the optimal control problem is the spectral factorization of 
an irrational transfer function. We highlight the basis for the spectral factorization algorithm 
together with associated computational issues arising in optimal regulator design. We also 
highlight options for implementation of wide band vibration control for flexible structures 
based on the sampled-data frequency response models. A simple flexible structure control 
example is considered to demonstrate the combined frequency response modeling and control 
algorithms. 

1 Introduction 

Frequency response methods offer several advantages for modeling the dynamics of small 
amplitude vibrations in flexible structures. Such models capture the distributed parameter 
dynamics of the elastic structural response as an irrational transfer function from localized 
actuation to localized deformation measurements. Interest in frequency response models can 
arise from a desire to predict modal frequencies with increased accuracy over that obtainable 
from finite element methods. The frequency domain approach is well suited to optimal control 
law synthesis with specific requirements for precision vibration supression and isolation. Most 
computational methods for optimal control synthesis available to design engineers focus on 
the manipulation of state space models. For flexible structure control, state space models 
are problematic since the question of model order required must be resolved as part of the 
optimal control computation. This paper reports progress in the development and testing of 
computational methods for design of precision control systems for mechanical structures with 

1 Work supported by SDIO and Air Force Wright Research and Development Center under contract F33615- 
88-C-3215. 


277 



C + 

open right half complex 5-plane, 9?e s > 0 

Hoo 

Hardy space of complex functions, analytic and essentially 
bounded in C + 

h 2 

Hardy space of complex functions, f(s ), analytic in C+and such 

that ; O \\fi s )\\ 2 d3 ) 1/2 < °o for « € C+. 

RH„ 

rational functions in 


Table 1.1: Notation 


elastic effects based on direct frequency response models. The frequency response models 
can arise from finite element analysis, transfer function methods, wave propagation models, 
and/or empirical measurements. Moreover, the computational approach offers a framework for 
integration of frequency response data from various modeling approachs which offer varying 
precision in different frequency bands. The current paper extends the efforts reported in [1]. 

We will use the following notation and conventions in this paper. The transpose of a 
column vector will be denoted as Tr X is the trace of the square matrix A", and j — 
yj— 1. A Laplace (resp. z) transform will normally be indicated by dependent variable; x(s) 
(resp. x(z)) 1 however, we often drop the explicit dependence where the meaning is clear 
from the context. The notation u*(,s) = u T (—s) will be frequently used, indicates 

the expectation of the random process x(f). In this work all random processes are assumed 
wide sense stationary and ergodic so that expectation can be replaced with ensemble average 
where convenient. The notation contained in Table 1.1 specifies the classes of transfer function 
models considered at various points. A rational function has a (partial fraction) expansion 
A(a) = {A(s)}+ + {A(s)}_ + {y4(5)}oo where {.}+ (resp. {-}-) is analytic in 9?es > 0-t.he 
causal part ($?es < 0-the anti-causal part) and {.}«> is the part associated with poles at 
infinity. Thus the operation {A(s)} + is causal projection of the frequency response model. 

In section 2 we provide an overview of frequency domain models and modern Wiener-Hopf 
design of multiloop control systems. We motivate the role of frequency response modeling and 
optimal control and identify critical computational steps required for the method. Section 3 
discusses a new approach to the required computations for Wiener-Hopf control which extend 
the algebraic constructions for rational transfer functions to certain irrational cases. We 
highlight the role of coprime factorization in design of distributed systems. Section 4 considers 
a simple, but nontrivial distributed parameter system design. Finally, in section 5 we discuss 
new options for real time control implementation suggested by the computational approach 
of section 3. 

2 Optimal Control of Frequency Response Models: Wiener-Hopf 
Design 

Frequency domain models have been used to articulate the full range of opportunities for 
feedback compensation for internal model stabilization. Algebraic constructions based on 
Laplace transform models of linear, time-invariant system dynamics have been used to describe 
alternatives for standard control computations and realizations for stabilizing controllers [2, 
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3]. This together with Wiener-Hopf optimization provides a general approach for resolving 
tradeoffs in regulator design where natural, frequency domain specifications for model-based, 
control performance are available. We remark that restricting attention to rational transfer 
function models in computational approachs to flexible structure control has been primarily 
motivated by convenience [1]. Specific results which extend the constructions of coprime 
factorization and internal stabilization to a certain class of irrational transfer functions have 
been obtained [4]. In the present effort we restrict attention to transfer functions in H^and 
meromorphic with the exception of a small number of right half plane poles. 


Optimal regulator design via Wiener-Hopf methods. Techniques for the solution of 
H 2 optimization problems in multiloop feedback systems have received considerable attention 
in the control theory literature for a number of years. A comprehensive approach to Wiener- 
Hopf design using transfer function models is given by Youla et al [5]. 

A general framework for resolution of tradeoffs in multiloop control design was recently 
outlined by Park and Bongiorno [6]. In general, control design involves the resolution of 
choices in the use of dynamic (feedback) compensation with respect to a nominal dynamic 
model of the system response to an n - vector control, u and an m- vector of exogenous system 
disturbances, e, as seen by p available sensors, y and t (possibly nonmeasurable) regulated 
variables. A frequency domain model for the control design problem (shown in Figure 2.1) 
can be expressed using Laplace transforms as, 



G(s) 


u(s)\ 

A 3 )) 


Gyu Gy t 

G zu G ze 



(2.1) 


The control architecture is assumed to involve feedback, 


u(s) = C(.s)y(4). (2.2) 

Then the closed loop compensation will alter the response of the system to disturbances as 
seen in terms of the regulated variables as 2 , 

z = [g zu cs,i] e 

Lr ze 

= G zxt C SG ye e -\~ G ze e (2.3) 

where the system closed loop sensitivity operator is S = [I + GyuC']” 1 . 

A major consideration in the classical methods of frequency domain design is closed loop 
stability. In such methods stability considerations must be continually evaluated (using root 
locus or Nyquist plots) as performance tradeoffs are evaluated. For single loop designs of rel- 
atively low order systems, classical frequency domain methods focus attention on the tradeoff 
between stability margins and performance. The modern approach is to use optimization to 

Suppressing dependence on the Laplace variable s. 


279 




Figure 2.1: General Multiloop Control System Design Problem 


resolve complicated engineering tradeoffs in multiloop design — subject to the constraint of 
system internal stability. 

The well known Youla parametrization of all stabilizing feedback controllers C which 
stabilize a given plant model was orginally derived for rational matrix transfer functions 
which have coprime factorizations over the ring of polynomials in the Laplace s- variable [5]. 
It is now understood that the construction goes through without modification for the ring of 
stable rational functions, RHa, which includes all rational transfer functions analytic in the 
closed right half plane including the point at infinity [7]. Thus under the assumption that the 
plant transfer function has right and left factorizations, 

Gyu = ND- 1 = Df'Nt, (2.4) 

coprime over RH^, then there exist A , Y, X(, Yt G RH,* such that 

D ( X t + N t Y t = 1, X D + YN = /, (2.5) 

and each controller, C, which obtains R = CS with R G RHo, can be parametrized by the 

factorization formulae, 

C = (X - KN t )~ l {Y + KD t ) (2.6) 

= (Yt + KD){ X t - KN )"* (2.7) 

for some K £ RH*,. The importance of this construction is that optimization procedures can 
be applied directly to the choice of K G RH* without concern for closed loop stability. This 
fact was first exploited by Youla et al[5] in the description of Wiener- Hopf optimal control for 
frequency response models. More recently, the parametrization has been utilized by Desoer 
and his students [8]. 

Our concern here is with H 2 optimization for a certain class of irrational transfer func- 
tions which are “pseudo- meromorphic” in the sense stable coprime factorizations exist; i.e., 
the constructions in (2.4)-(2.7) obtain closed loop stability with N, D, N ( , D(,X, Y, Xt, Yt,K G 
Hoc[l, 4]. In the research reported herein we avoid algebraic constructions related to stable, 
coprime factorization of irrational transfer functions as considered by Desoer [7] and instead, 
focus on the development of numerical algor thims for approximating the frequency response 
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of the required objects. Such transfer functions can be adequately approximated over fi- 
nite frequency ranges by (rather high order) rational transfer functions. Models of this type 
arise in the study of vibrations of multibody systems with flexible interactions [9] and wave 
propagation in flexible mechanical structures [10, 11]. 

PSD modeling of control processes for performance specification. A wide class of 
standard control design problems including symultaneous requirements for tracking, distur- 
bance rejection and accomodation, etc. can be represented in the form of the general linear, 
time-invariant regulator problem of Fig. 2.1 and (2.1)-(2.2) where the objective is to choose 
a controller which stabilizes the closed loop system and minimizes a performance criterion in 
the form, 

J = rr[Q(.)fi(.)] is, (2.8) 

where P z is an effective Power Spectral Density (PSD) of the regulated variables, 2 . 

The representation of regulation performance in terms of PSD is quite practical for a va- 
riety of design problems arising in multiloop systems and provides a frequency dependent 
specification of control performance consistent with design requirements for vibration rejec- 
tion. One can extend the significance of PSD modeling to include a wide range of practical 
design considerations. The regulation PSD, P z , can be related to modeling assumptions on 
the exogenous inputs in terms of the closed loop transfer functions; 

P z = [ G UZ R , /] ^ P e (Gye,, <?«.)• (2 9) 

PSD models of exogenous inputs may include deterministic transient effects together with 
steady state stochastic PSD, $ e ; viz., 

P'(s) = A!£;{e(s)e.(s)} + A 2 #«(a), (210) 

Formulation of the performance objective, J, may include real, positive, A;, i = 1,2 which 
permit scaling relative importance of steady state and transient considerations to the com- 
posite performance and C?(s) is included to permit frequency weighting. PSD modeling has 
recently received increased emphasis in the study of vibration control in acoustic regimes [11]. 

Park and Bongiorno [6] also highlight the use of PSD models for minimizing closed loop 
system sensitivity to model uncertainty. Let the system model uncertainty be given as a 
frequency dependent, additive perturbation, G := (7+A, which can be expressed in partitioned 
form as, 



Following [6] an effective model uncertainty PSD, Pa = 2?{AA„}, can be reflected to the 
system regulated outputs, and via superposition, a composite performance objective of the 
form, 

j =^C Tr { Q ^- R ^[tz S:D* 
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is obtained for optimal regulator design, where the system model uncertainty can be obtained 
commensurate with the performance specifications as an effective disturbance PSD, 


<J> 

*yy ^v z 
$ d> 

* Zli T zz 



P e (G ye .,G t '*) + nP a - 


( 2 . 11 ) 


Here the partitioned terms can be expressed in terms of a priori modeling assumptions; [6] 

= G ye P e G y '.+n(E{A vu (&ME{A ye (& y ,).}), (2.12) 

= G ye P e G ze * + p(£J{Ay U (A zu ),}£/{A ye (A ze ),}) = (213) 

= G ze P e G„,+ti(E{A zu (A zu ),}E{A ze (A ze ).}). (2.14) 


It is by now widely recognized that frequency domain response considerations are ex- 
tremely important for robust control design and that performance objectives formulated in 
the frequency domain are important tools for resolving design tradeoffs of relevance to prac- 
tical design problems. However, the common wisdom is that state space modeling offers the 
most reliable numerical framework for the computational problems which arise in optimal 
regulator design. The Wiener-Hopf approach identifies the solution for the optimal controller 
in an explicit form which highlights the role of the algebraic constructions generic to stabiliza- 
tion and the quantitative computations required for identifying an optimal controller. Thus 
given the system architecture (2.1)— (2.2), appropriately chosen stable coprime factors for the 
plant (2.4), a nominal stabilizing controller given in terms of its coprime factors as solutions of 
the Diophantine relations (2.5), and performance PSD’s (2.12)-(2.14), then an optimal closed 
loop system response is obtained (assuming a solution exists) by the formula, 

R = DA- 1 ({A D-'YU}- - {A: 1 DMtu.QQzyDtM: 1 }^ fT 1 #,. (2.15) 

The explicit form given here depends on operations of causal projection and the solution of 
two causal, spectral factorizations; 

D.G ZU .QG ZU D = A.A (2.16) 

D($yyDl* = (2.17) 

with A, A -1 , S7, Q -1 € H,*,. The required controller can then be obtained in the explicit form, 
C = (I - RGyv)- 1 R. 

The computational steps required to identify candidate optimal control solutions for the 
regulator problem include: 1) stable coprime factorization (as in (2.4), 2) identify candidate 
solution to Diophantine relations (2.5), 3) causal spectral factorization, and 4) causal pro- 
jection. We contend that such computations can be effectively supported (in finite precision 
arithmetic) by obtaining state space realizations [3] only for relatively low order , rational trans- 
fer functions. In the sequel, we specifically avoid such an approach since we are ultimately 
concerned with the approximate solution of large (or even infinite) dimensional models. 


3 Frequency Response Computations for Optimal Control. 

Our approach to optimal control computation is motivated by distributed parameter models 
which arise in flexible structure control. The approach we have in mind is based on sampling 
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and interpolation of the frequency response models for the system. The choice of sampling and 
the resulting high order, rational approximations are obtained in the context of the optimal 
control problem as summarized above. 


A computational approach to spectral factorization. Recall that a transfer function 
H(s) £ H 2 fl Hqo has a unique spectral factorization H(s) — F($)F*(s) with F £ if: 

1. H(s ) = H(s ); i.e. , H(s) is the transform of a real- valued function h(t). 

2. H(s) ~ i.e., H(s) is “para-hermittian” . 

3. H(s) is of normal rank; i.e., full rank almost everywhere in C. 

4. H(iu>) is positive, semi-definite and bounded for u; £ R. 

To see that causal projection is a closely related problem consider the following. If H(s) is 
scalar, then with $(s) = In H(s) we obtain 

*W = {*(.)}.+{*(.)}-, (3.1) 

= lnF(») + lnF.( S ), (3.2) 


so that the causal, spectral factorization is related to causal projection via the logarithmic 
transformation; F(s) — exp{\n H(s)}+. 

Our goal is to obtain numerically stable approximations to these related problems for 
transfer functions in H^. For application to precision control of flexible structures we require 
wide band frequency domain models so that even rational approximations will be of relatively 
high order. An approach to model order reduction which has recently received attention in the 
literature is based on Fourier series approximation of irrational frequency responses [12] in Hoq. 
Our approach to computations for such models is also based on sampling and interpolation of 
the spectrum, but is motivated by computational requirements for Wiener-Hopf optimization. 
From the above discussion of causal projection we motivate a class of algorithms of interest 
from basic properties of the Hilbert transforms applied to the frequency response Recall 

that the Hilbert transform of a time signal f(t) is defined as a convolution; f(t) = dr 

and it’s Fourier transform has the property, 


/(«) 



j/M, 


(X) > 0 

o> > 0 


The inverse Fourier transform of $ is — jsgn(t)<f>(t) where <f>{t ) is the inverse Fourier transform 
of $(a;). A consequence is that the casual projection can be obtained as 

{♦(«)}+ = + 

In previous studies we reported computational algorithms for causal projection and scalar 
spectral factorization by numerical evaluation of the Hilbert transform integral. Computa- 
tional cost was high due to the fact that the Hilbert transform integral is convergent only in 
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the Cauchy principal value sense [13]. An alternate method for causal projection and spectral 
factorization was considered in [14] based on sampling and interpolation of the system fre- 
quency response. The algorithm developed in [14] employs results of Stenger [15] on numerical 
solution of Wiener-Hopf integrals by sampling and interpolation. Details of the algorithm used 
for the current studies and computer implementation are given in [14]. 

In the multiloop, optimal regulator design problem we require the solution of two matrix 
spectral factorization problems (analogous to the solution of control and filtering Riccati 
equations for time domain models). The computational approach exploited in the current 
study is based on a Newton-Raphson iteration for the matrix causal spectral factor; 

F n+1 (iu,) := {[F n *(^)]- 1 i/(iu;)[F n (^r 1 } + F n (*w). (3.3) 

The recursion (3.3) can be replaced with a numerically well conditioned problem by iteration 
on the inverse spectral factor; 

== w (/ + {ra- 1 h [F„r i - /} + )'* <3.4) 

By initializing with P 0 (an m x m diagonal matrix) with diagonal elements equal to the spectral 
factors of the diagonal elements of H the second term of (3.4) remains a perturbation of the 
identity (since [F*]~ l HlFn}" 1 — I — > 0) which regularizes the computations. The algorithm 
used in this work is based on that reported in [14] and is a modified form of the method 
reported in [16]. 

Computation of stable coprime factorizations for flexible structure models. Simple 
models of structural components with elastic effects typically lead to transfer functions in 
HooOnce realistic damping models are included. Linear vibration models of more complex 
structures arising in aerospace applications usually will have transfer function models with 
only a finite number of poles in the closed right half plane. Restricting attention to such 
transfer functions we indicate a simple procedure for coprime factorization over 

Let S C Hqo be a set of transfer functions analytic in a half plane including C+. Under 
the above assumption any such transfer function P(s) can be expressed in the form, 

P{s) = Ps{s) + Ps(3) (3.5) 

where 6 <5 C and Pg is rational and analytic in the complement of S with (a finite 
number of) poles outside S. A stable coprime factorization can be readily obtained for the 
(typically low order) transfer function as, Pg = N r D~ l , by well known state space construc- 
tions [3]. Then P has stable coprime factorization, 

P = N r D = [N r - P s Dr]b;\ (3.6) 

where N r , D r are 5-stable. The separation of terms in (3.5) is readily carried out given P(s) 
by computing the residues of the finite number of unstable poles contributing to Pg. 
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4 Control Computations for an Elastic Structure 


To illustrate the computational approach for a simple elastic structure we consider the simply 
supported Euler beam with torque control at one end. The beam lateral deformation is given 
by y(t,z) with 0 < z < L and has dynamics described by the dimensionless PDE; 


dy 2 

8t 2 


~2C 


<9 3 y 

dtdz 2 


d A y 


(4.1) 


with boundary conditions at z = 0, 

y(t, o) = o, 

and at z — L, 


cPy 

dz* 


= o, 


2=0 




= r. 


i=L 


with the control moment, r, applied at the right hand end of the beam. The transfer function 
(r — > y) for beam control is 


G yu (s,z) = L 2 


sin A x sinh A 2 f — sin A x jr sinh A 2 
(Af + Al)sin(A,)sinh(A 2 ) ’ 


(4.2) 


where A, = (— ( + iy/1 — C 2 ) s ^ 2 > A| = (( -f iy / 1 — ( ^)sL 2 , L is the beam length, ( is the 
damping factor, and z is the observation point to be regulated on the beam. The transfer 
function is meromorphic, and G vu (s, z ) & Hoc for any 0 < z < L. 

The regulator problem considered arises from a requirement for asymptotic rejection of 
constant load disturbances at a point £ = 0.7. For the current numerical studies we take 
L = 10., and the effective damping ratio, ( = 0.01. Stable coprime factorization is trivial 
and we take N r = N( = G yu , D r = D( = 1. Exogenous inputs here include the output 
load disturbance d and measurement noise model n and are described by their effective PSD 
models representing constant (step) load disturbance and narrowband sensor noise as shown in 
Figure 4.2. The frequency response of Gy^ is shown in Figure 4.1 with 1024 uniform frequency 
samples over a bandwidth of 0 < u) < 100. Clearly, the frequency response is irrational and 
no obvious rational approximation is evident. 

The optimal control design is regulation of the beam deflection at zjL = .7 and the 
performance objective is given as, 

J = : / Tr{$ u + /t$ u } ds 

Z7TJ J - joo 

where the tracking cost is modeled by PSD, $ y , and the control saturation PSD is $ u = 
£{tiu„}. Then given a constraint on the control power the scalar fi > 0 plays the role of a 
Lagrange multiplier for the optimal design. In this case the required spectral factors; 


(G yU* Gyu + ft) — A„A 

{Gyd^dPyd- + $n) = DD, 


(4.3) 

(4.4) 
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Figure 4.1: Frequency response for pinned- pinned beam control. 

were obtained for p = .1 by the frequency sampled procedure and are displayed in Figure 4.3. 
We remark that the computations of the indicated spectral factors effectively replace the com- 
putational step of solving a pair Riccati matrix equation for the control (resp. filter) problem 
typically encountered in state-space methods for control design. For distributed parameter 
systems, solution of the Riccati equation (a PDE) requires discretization which is accom- 
plished using the current algorithms by sampling and interpolation of the frequency response. 
Thus numerical precision is concentrated over frequency bands significant for the given control 
problem and with sampling under direct control of the designer. The solution obtained is ef- 
fectively a high order rational approximation of the optimal solution with frequency response 
interpolation points chosen by the design engineer. The optimal controller frequency response 
thus obtained is shown for p = .1 in Figure 4.4. 

5 Frequency Sampling Filters for Real Time Control Implemen- 
tation 

The frequency response computations for Wiener-Hopf control outlined and illustrated in the 
previous sections identify various frequency sampled approximations to the ideal, possible 
irrational frequency response for the desired optimal controller. Bandwidth and sampling can 
be chosen by the design engineer to represent specific concerns based on models and/or control 
performance. The frequency sampled computations obtain a specification for the frequency 
response of the ideal (optimal) controller via its sampled representation. The design engineer 
now has several options for implementing the controller depending on available hardware. In 
contrast to the state space approach for finite dimensional systems, several new realization 
opportunities are suggested by the frequency sampling approach. 
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Figure 4.2: PSD for disturbance and sensor noise inputs for beam control. 


A principal concern in implementation of high precision control laws for flexible structure 
control is the order of the realization required for the online controller. The controller or- 
der is usually taken to mean the dimension of the state variable realization of the transfer 
function C'(s) which will be implemented for realtime control. Implementation using analog 
components of high order models is limited by circuit complexity, reliability and cost. As 
a result considerable effort has been expended in methods for model order reduction. One 
approach to controller realization which follows from the frequency sampled computations of 
the previous section is to compute reduced order, continuous time, state space realizations for 
the controller by techniques such as in [12]. 

Digital computer implementations are primarily limited by computational speed and al- 
gorithm complexity effecting the ultimate obtainable sampling rate and considerations for 
reduced order realization of the controller may be required. However, the emergence of special- 
ized computer hardware implemented in VLSI single chip circuits for digital signal processing 
opens new opportunities for realization of realtime control for flexible structures. We prefer 
to consider realization options for the optimal controller in discrete time for implemenation 
on a digital computer. Realization of the controller specified by its frequency samples can be 
obtained using a FIR digital filter implementation. 

Given the specified frequency samples obtained for the optimal controller, 

c k = C(j u/fc), 

at frequencies, a > k = kwBw/N, where ujbw is bandwidth and N the number of uniformly 
spaced frequency samples we describe the digital filter realization using z-transforms. With 
discrete time sampling rate u>, > 2 u>bw the frequency samples correspond to interpolation 
points in the z-plane given by 3 , z h = e jk2ir/Nl , for k = 0, . . . , N — 1. The z transform which 

3 Bandwidth and sampling requirements would typically require padding the sequence of frequency samples 
of length N with Ni — N zero values to avoid aliasing. 
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Figure 4.3: Spectral Factors for pinned-pinned beam control. 
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realizes the frequency sampling filter is 


N - 1 


C(z) = £ C,*i(z) 

fc =0 


where the interpolating functions are, 

1 - z-" 

Fk[2) ~ N{ 1 - e^r>)’ 

with & = 0, . . . , N — 1. A standard computation shows that the frequency sampling filter has 
transfer function 

i _ r -N N - 1 /-< N-\ 

= = (*•») 

where the coefficients, 

(5.2) 

JV fc =0 

for t = 0, ...,N — 1, are the Inverse Discrete Fourier Transform (IDFT) of the sequence 
Co, . . . , Cjv_i. The final form in (5.1) shows that the realization is a FIR realization. Such 
realizations are nonrecursive and are efficiently implemented using high speed, single chip DSP 
processors which utilize highly pipelined architectures to achieve high throughput. 

6 Conclusions and Directions 

Wiener- Hopf optimization of frequency domain models has been shown to offer significant 
advantages for computation of precision controllers for irrational transfer functions arising 
in control of flexible structures. Computational algorithms for causal spectral factorization 
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and causal projection can be implemented based on frequency sampled representation of the 
model response. Such models can be obtained from transfer function models or from frequency 
response measurements of the controlled structure. Computations based on frequency response 
sampling have been demonstrated for irrational transfer function models arising in the control 
of flexible structures. 

Requirements for precision control will involve frequency response models which are char- 
acterized by a large number of flexible modes within the control bandwidth. However, for 
control of relatively large, flexible space structures control bandwidth and resulting sampling 
requirements for discrete time control implementations are well within the state-of-the-art for 
high speed digital computers. Frequency sampling filters based on nonrecursive implementa- 
tions can be efficiently implemented in modern DSP single chip processors for realtime control 
of such systems. 

Application of FIR realizations for realtime, closed loop control have not received much 
consideration in the literature primarily due to increased phase lag by comparison with a 
recursive realization. However, the rapidly developing technology for realtime DSP using 
special purpose architectures offers throughput capabilities which may reduce the achievable 
computational delay to within acceptable limits for certain applications. In such cases, high 
order realizations may be feasible using nonrecursive impleinenations which cannot be realized 
by recursive methods. 
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